Extinction of quasiparticle interference in underdoped cuprates with coexisting order 
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Recent scanning tunnelling spectroscopy measurements [Y. Koksaka et a/., Nature 454, 1072 
(2008)] have shown that dispersing quasiparticle interference peaks in Fourier transformed conduc- 
tance maps disappear as the bias voltage exceeds a certain threshold corresponding to the coincidence 
of the contour of constant quasiparticle energy with the antiferromagnetic zone boundary. Here we 
argue that this is caused by quasistatic short-range coexisting order present in the d-wave super- 
conducting phase, and that the most likely origin of this order is disorder-induced incommensurate 
antiferromagnetism. We show explicitly how the peaks are extinguished in the related situation 
with coexisting long-range antiferromagnetic order, and discuss the connection with the realistic 
disordered case. Since it is the localized quasiparticle interference peaks rather than the underlying 
antinodal states themselves which are destroyed at a critical bias, our proposal resolves a conflict 
between scanning tunneling spectroscopy and photoemission regarding the nature of these states. 

PACS numbers: 74.20.-z, 74.25.Jb, 74.50.+r, 74.72.-h 



I. INTRODUCTION 

The understanding of how a Mott insulator with local- 
ized states becomes a metal as one gradually increases the 
carrier concentration remains one of the main challenges 
of condensed matter physics. This question may be inti- 
mately connected with the so-called nodal- antinodal di- 
chotomy (sharp quasiparticles in the nodal region versus 
broad, gapped antinodal "quasiparticles") observed by 
angular resolved photoemission spectroscopy (ARPES) 
on underdoped cuprate materials. ^'^ Tunnelling spec- 
troscopy has also been used to probe states in different 
regions of momentum space by the help of Fourier trans- 
form scanning tunneling spectroscopy (FT-STS).^~^ In 
particular, it has been argued that the apparent incoher- 
ent antinodal states have their origin in the emergence of 
charge ordered regions in the underdoped regime.^ 

In the d-wave superconducting (dSC) state near opti- 
mal doping, quasiparticle interference (QPI) observed by 
FT-STS is dominated by peaks at well-defined wavevec- 
tors which agree well with those predicted by the 
so-called octet model. ^'^'^ The contours of constant en- 
ergy (CCE) of the Bogoliubov quasiparticle dispersion 
are generally shaped like curved ellipses ("bananas") 
centered at the nodal points; seven different nonzero 
wavevectors connect the tips of these bananas. The 
dispersion of the peaks at q^ allows one to extract the 
shape of the underlying Fermi surface^^ as well as the 
momentum dependence of the superconducting energy 
gap.^ A quantitative understanding of the amplitude and 
width of the peaks, however, is not straightforward to ob- 
tain, and depends rather sensitively on the nature of the 
scattering medium. A point-like scatterer, for example, 
in an otherwise homogeneous dSC leads to a landscape of 
some spot-like and some arc-like dispersive features close 
to qi in the FT-STS images^' "'^"'^ "'^^ whereas experimental 
data appear mostly spot-like. Furthermore, the weights 



of the various peaks calculated in simple theories differ 
from experiment. There have been several theoretical at- 
tempts to remedy this situation by including more realis- 
tic models for the disorder present in Bi2Sr2CaCu2 08+5 
(BSCCO). In particular, as realized recently, it appears 
important to include the gap inhomogeneity arising from 
the dopant atoms. ^^"^^ Models have also included the ex- 
tended Coulomb potential arising in this material from 
partly screened Bi^Sr substitutional disorder and the 
oxygen dopant atoms. Nevertheless, a complete quan- 
titative description of the FT-STS patterns is still lacking 
in the superconducting state. 

FT-STS has also been used to probe the pseudogap 
state in the underdoped regime where non- dispersive 
(bias-independent) quasi-periodic conductance modula- 
tions has been identified both in BSCCO^'^^'^^ and 
Ca2-,,Na^Cu02Cl2 (Na-CCOC).^^ The origin of these 
peaks remains unknown at present but may be caused 
by short-range charge order, possibly connected to the 
existence of nested segments of the Fermi surface near 
the antinodal regions of the Brillouin zone.^^ In BSCCO 
near optimal doping non-dispersive peaks have also been 
discussed in terms of pinned and disorder-induced charge 
order. ^^"^^ In Na-CCOC it has been proposed that 
phonons play a crucial role in stabilizing a d-wdive charge 
density wave order, or a surface transition to a com- 
mensurate charge density wave state.^^ At present, how- 
ever, it remains controversial whether true charge order- 
ing is required for describing the non-dispersive LDOS 
modulations . ^^^^ 

Recently, new developments in the FT-STS technique 
allowed for further detailed exploration of the electronic 
properties of underdoped Na-CCOC and BSCCO. For 
example, it was argued that tip-elevation errors can be 
avoided by studying the conductance-ratio Z{r^E = 
eV) = ^(r, V)/g{T^ —V) where V is the bias voltage and g 
the conductance, and the detailed properties of the LDOS 
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modulations were investigated in this regime as well.^^^^ 
It was found that irrespective of the doping level, an "ex- 
tinction line" exists in momentum space, beyond which 
most of the dispersing FT-STS peaks (q2,q3,q6 5q7) dis- 
appear, to be replaced by a reduced set (qi,q5) of roughly 
non-dispersive peaks. This extinction line is doping in- 
dependent and coincides at all doping levels with the 
antiferromagnetic (AF) zone boundary [lines joining the 
points (0, ±7r) and (±7r,0)]. At energies below the scale 
Ao where the CCE first touch the AF zone boundary, 
the FT-STS response is similar to the dispersing Bogoli- 
ubov quasiparticles obeying the octet model. At energies 
above Aq, the response becomes highly spatially inho- 
mogeneous and appears dominated by pseudogap excita- 
tions. 

One popular picture of the nodal-antinodal dichotomy 
invokes intense scattering with momentum transfer near 
(7r,7r) which broadens states near the antinodal points. 
The problem with this picture in the superconducting 
state, however, is that the phase space for scattering is 
smallest at precisely these points of momentum space 
because the (i-wave gap is largest there. Thus Graser 
et al.^^ calculated the spin fluctuation spectrum within 
an RPA-type formalism and used it to determine the 
lifetime of states near the node and antinode of a dSC 
phase. This same framework produces a good description 
of the resonant magnetic response near (tt, tt) as mea- 
sured by neutron scattering. Nevertheless their results, 
which were consistent with earlier work on quasiparticle 
lifetimes, ^^^^ imply that inelastic scattering of the con- 
ventional itinerant spin-fluctuation type cannot severely 
broaden quasiparticle states in the superconducting state 
with momenta near the antinodes. This is of course con- 
sistent with ARPES experiments, which find broad but 

well-defined antinodal peaks in the superconducting state 
of BSCCO.i'43,44 

One aspect of the physics of some underdoped cuprate 
materials which is left out of the conventional spin- 
fluctuation scattering analysis is the possibility of addi- 
tional order coexisting with the dSC state at low tem- 
peratures. Static stripe order has been observed in sev- 
eral high- Tc materials and appears especially pronounced 
near 1/8 doping;^^ in the La2-a:Baa^Cu04 (LEGO) sys- 
tem, for example, charge order appears around 50K and 
persists to lower temperatures, where it coexists with 
spin order. In addition, /iSR has consistently reported 
so-called "cluster spin glass" (GSG) signatures of frozen 
magnetic order all over the underdoped cuprate phase 
diagram of La2-:,Sr^Gu04 (LSGO) and BSGGO,'^^''^^ 
generally attributed to disorder present in signficant 
amounts due to the way in which these materials are 
doped. In both LSGO and BSGGO, the nodal-antinodal 
dichotomy is also observed in ARPES measurements.^^ 
These observations suggest that quasiparticle scattering 
from short-range coexisting order might play an impor- 
tant role in explaining the extinction of the QPI peaks 
in the experiment by Kohsaka et al.^^ 

There is no complete consensus on the origin of these 



ordering phenomena. One general notion is that disorder 
can pin fluctuating order while still reflecting the intrin- 
sic correlations of the pure system. Whether the level of 
disorder present in these intrinsically disordered systems 
is too large to justify such an assumption is not clear. 
However, various concrete models of pinned fluctuating 
stripes have been proposed^ which resemble experi- 
ment in qualitative ways. An alternative starting point 
to understand the GSG phase assumes that dopants nu- 
cleate droplets of staggered order, which then interfere 
constructively to create quasi- long- range order. ^^'^^ 

These studies are related to the important question 
of whether AF correlations coexisting with preformed 
Gooper pairs are adequate to explain the pseudogap 
phase, or whether other types of ordering phenomena 
occur instead, or in addition. ^^'^''"^^ Models of a dis- 
ordered AF phase coexisting with c?-wave pairs can re- 
produce a number of known experimental results: the 
magnetic correlations have been shown to protect the 
nodal quasiparticles;^^'^^ reproduce the Fermi arcs and 
nodal-antinodal dichotomy;^^ as well as the tempera- 
ture dependence of the superfluid density^^ and thermal 
conductivity.^^'^^ 

Although the coexisting order in these systems is short- 
range, we know from neutron measurements in under- 
doped LSGO^^ that correlation lengths can be quite long, 
of order 100 lattice spacings. It may therefore not be a 
bad starting point for the study of the effect of these cor- 
relations on QPI to assume a long-range ordered state. 
Here, we propose a concrete model for this effect by as- 
suming the presence of a long-range AF state coexist- 
ing with superconductivity, and then investigating the 
consequences for the LDOS modulations. Even though 
the underdoped materials which exhibit magnetic order- 
ing are characterized by incommensurate antiferromag- 
netism [with incommensurate peaks away from (7r,7r)], 
we assume (7r,7r) ordering for simplicity. We focus on 
the origin of the extinction line, and not the high-energy 
LDOS modulations which require more realistic disorder 
models, possibly including charge ordering. In order to 
clearly elucidate the role of the AF order we consider 
for simplicity a single point-like scatterer. The results 
presented below should be important for understanding 
future FT-STS modelling using more realistic disorder 
configurations. 



II. FORMALISM 

The FT-STS signal of a disordered dSG has been dis- 
cussed rather extensively by theoretical mo dels. ^'-^-^"-^^ In 
the present case of coexisting AF order, the translational 
symmetry is broken, and the formalism is very similar to 
the (i-density wave approach discussed in e.g. Ref. 65. 
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The Hamiltonian reads 
H = ^[e(k)-M]4,Cko 



(1) 



+ H.C., 



where cj^^ creates an electron with momentum k and 
spin (7, M is the magnetization and Q = (tt, tt) is the 
AF ordering vector. We work in units where the lattice 
constant a = 1. The superconducting c^-wave gap func- 
tion is A(k) = A {coskx — cos ky) /2, and the quasipar- 
ticle dispersion is e(k) = ei(k) + 62 (k), where ei(k) = 
—2t (cos kx + cos ky) and 62 (k) = —At' cos kx cos ky — 



2t" {cos2kx -\- cos2ky). It is convenient to split up 
the normal state band in this form due to the differ- 
ent symmetry properties of ei(k) and 62 (k) with re- 
spect to momentum shifts of the AF ordering vector 
Q; ei(k + Q) = -ei(k) and e2(k + Q) = e2(k). In 
terms of the following generalized Nambu spinor ip^. = 

{cj^p cj^^Q|, c_kX5 c_k-Q|}, we can write the Hamilto- 
nian in the form 



(2) 



where the sum is restricted to the reduced Brillouin zone 
(RBZ), \kx \ + \ky\ < TT, and A(k) is given by 



ei(k) + e2(k) -/X M A(k) \ 

M -ei(k) + e2(k)-M -A(k) 

A*(k) -ei(k)-e2(k)+M M 

-A*{k) M ei{k) -e2{k)+/i/ 



(3) 



The eigenvalues zbE'i^2(k) of A{k) are given by 



^l,2(k) = y (^[62(k) - /i] ± y^6Kk)+ M2j + A2(k), 

(4) 

which for A(k) = reduces to E^(k) = [e2(k) - 
//] ± ^e?(k) + and for M = reduces to ^(k) = 
±v'[ei(k) + e2(k) - + A2(k). The Greens function 
Go(k, iojn) of the pure system is obtained from the equa- 
tion 

Go(k, icJny^ = iujnl - A{k), (5) 

where / denotes the 4x4 identity matrix. 
In the presence of an impurity term 

Himp= 4^(k,k0^k^ (6) 

k,k'Gi?BZ 

the impurity-contribution to the full Greens function 
G(k, k', ic^n) is given by 

G(k, k^ iujn) = Go(k, icjn)^(k, k^ i6Jn)Go(k^ icj^), (7) 

where 

^(k,k^^u;0 = y(k,kO+ (8) 
J2 V{k, k'OGo(k^ iuJn)T{k'\ k^ iujn). 

For a point-like impurity V{k^k') [and T(k, k', io;^)] be- 
come independent of k and k^ Specifically, a nonmag- 



netic (5-function scatterer takes the form 

/ 1 1 \ 

y(k,k') = y hj_o, (9) 

\0 -1 -1/ 

The change in the LDOS from the pure phase 
^A/'(q, cj), is given by^^ 

^iV(q,w) = ^ E 5(k,q,a;), (10) 

\<ieRBZ 

where ^(k, q, cc;) is defined as follows. Let k' = k + q. If 
k' is in the RBZ, then 

4 

5(k, q, io) = E[Gii(k, k', Siio) - G*,(k', k, s^uj)], (11) 

i=l 

where Si = 1 for for the particle- hole sector i = 1, 2 and 
Si = —1 for the hole-particle sector i = 3,4. If k^s not 
in the RBZ, let k'' = k + q - Q. For this case 

^(k, q, u;) = ^ [Gi,i+i (k, k'', siu) - G* (k'', k, siu) 

2 = 1,3 

+ G,+i,,(k,k^8,c^) -G*+i,,(k^k,5,c^)].(12) 

Here, G(k, k',c<;) is obtained by usual analytical contin- 
uation iujn ^ (jO -\- iO+ of G(k, k', zcj^). Below we in- 
troduce a finite lifetime broadening r] = 0.02t such that 
iuon ^ (jO -\- if], and the summation over the RBZ is per- 
formed using a 600 x 600 mesh. 

We use a band dispersion with t = 1.0, = — 0.3,t'' = 
0.1 and /i = —1.25, which yields the normal state Fermi 
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FIG. 1: (Color online) (a) Contours of constant quasiparticle 
energy for the pure d-wave superconductor with A = 0.6t at 
energies cv/t = 0.0,0.075,0.225,0.375,0.5,0.57. Also shown 
are the underlying normal state Fermi surface, the AF zone 
boundary, and the seven distinct nonzero wavevectors con- 
necting the banana tips. 



surface shown in Fig. 1. In addition, when studying the 
superconducting state we set A = 0.6t. In the follow- 
ing we focus the discussion on energies below this gap 
energy, which is taken unrealistically large for numerical 
purposes in order to highlight features at low energies. 
The impurity potential V = 0.1 1 is chosen to be weak in 
the sense that it does not produce any low-energy reso- 
nant states. For simplicity, we ignore any spatial struc- 
ture of the local Wannier orbitals, rendering all results 
periodic in momentum space with respect to reciprocal 
lattice vectors. 
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FIG. 2: (Color online) QPI maps versus Qx/tt and Qy/Tr in the 
dSC phase in the presence of single point-like potential scat- 
terer of strength V — O.lt. Energies uj are as shown in the 
Figure titles. Some dispersing spot-like features correspond- 
ing to octet vectors are easy to identify. For example, we 
have circled the qr peak [see Fig. 1]. 



III. SUPERCONDUCTING PHASE 



The QPI patterns in the pure dSC phase have been 
thoroughly discussed in the literature, and we will not 
dwell on them here. However, in order to discuss the 
effect of AF on the QPI, we show briefly some results 
for the pure dSC phase in this section. Figure 1 displays 
some typical CCE exhibiting the usual banana shaped 
form, and centered at the nodal points (±0.36, ±0.36)7r. 
The wavevectors connecting the tips of these bananas, 
reveal where peaks in the FT-STS maps are expected, 
although matrix elements consisting of certain combina- 
tions of coherence factors are important for this simple 
picture to hold.^'^^ Here we calculate the Fourier trans- 
form density of states \SN\ with SN given by Eq.(lO), 
and refer to it as a QPI map. Figure 2 shows QPI maps 



versus Qx and Qy at representative fixed energies inside 
the gap. For clarity, we have circled the qr peak which 
is positioned along the (110) direction and disperses to 
higher momenta with increasing energy as expected from 
Fig. 1. Line cuts along the nodal and antinodal direc- 
tions for the pure dSC phase are shown in Fig. 3 for 
both positive and negative energies, exhibiting clearly 
the dispersive QPI peaks discussed in the literature. The 
dispersion agrees well with the octet model, ^ which as- 
sumes that the scattering peaks are determined entirely 
by the wavevectors connecting each banana tip with all 
the others. It is easy to identify in Fig. 3 the octet peaks 
qa and qr in the (110) cut, and the qi and qs peaks in 
the (100) cut. 
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FIG. 3: QPI line cuts along the (110) (a,b) and (100) 
(c,d) directions for the pure dSC phase. Panels (a,c) 
[(b,d)] correspond to positive [negative] energies at |cj|/t = 
0.01,0.075,0.15,0.225,0.3,0.375,0.45,0.5 (bottom to top). 
For clarity the curves are displaced by 0.05. 



IV. ANTIFERROMAGNETIC PHASE 



In the pure dSC state, the spectral weight at the Fermi 
energy A(k, = 0) consists of four nodal points, which 
may be smeared slightly by disorder. On the other hand, 
in the pure metallic AF case, the spectral function can 
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FIG. 4: (Color online) Spectral weight A(k,oo — 0) in the 
pure AF phase with M = 0.5t (a), and M = t (b). The lower 
two panels (c,d) show CCE at the same energies as in Fig. 1 
for M 0.5t (c), and M = t (d). 



FIG. 5: (Color online) Same as Fig. 2 but for the pure AF 
phase with M — t. 



be written as 

A{k,co) =v\k)S[co-E]^{k)]W{mco-El,{k)i (13) 



with u^{k) = 1/2[1 + ei(k)/v'M2 + e2(k)] and ^^(k) = 
l/2[l-ei(k)/v'M2 + e^(k)]. In Fig. 4 (a,b) we show the 
spectral weight A(k, c<; = 0) in a case with M = 0.5t 
(a) and M = t {h). In both cases only the E]^ band 
crosses the Fermi level. As seen in Fig. 4(a,b), the outer 
ring of the Fermi pocket is washed out because of the 
'^^(k) coherence factor. This is caused by the unit cell 
doubling in the AF state, and a similar effect happens 
in e.g. the d-density wave scenario. The checkerboard 
charge-order scenario for the pseudogap phase also repro- 
duces a Fermi arc due to the difference in the coherence 
factors between the inner and outer parts of the arc.^^'^*" 
Disordered antiferromagnetism will further enhance this 
apparent spectral weight suppression on the outside of 
the pockets. 

For the case with M = twe show in Fig. 5 typical QPI 
maps in the pure AF phase. The QPI images are dom- 
inated by arcs of scattering intensity, rather than spots. 
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The reason why distinct, isolated scattering wavevectors 
do not occur in the pure AF case is caused be the differ- 
ent coherence factors in this case. As discussed in Ref. 12 
in the case of intra-nodal scattering, the origin of peaks 
in the QPI maps in the pure dSC phase is caused by the 
dSC coherence factors which conspire to significantly en- 
hance the weight near the tips of the CCE resulting is 
a significant enhancement of the QPI response localized 
near qy. In the pure AF phase, on the other hand, the 
different coherence factors cause the weight to be evenly 
distributed along the CCE resulting in arc-like character- 
istic QPI features. We note that the momentum-resolved 
density of states, which is qualitatively similar (at least 
at negative energies) in the two phases, is not the cause 
of this qualitative difference in the QPI maps. 

In Fig. 6 we show line cuts along the nodal and antin- 
odal directions at both positive and negative energies. 
We can understand the cusps in these line cuts from Fig. 
4(d): consider for example the results presented in Fig. 
6(a) showing cuts at positive energies along the (110) di- 
rection. The CCE clearly form ellipses centered around 
each nodal region as seen from Fig. 4(d). Scattering can 
be inter-ellipse or intra-ellipse, and leads to four cusps in 
agreement with Fig. 6(a). Since at positive energies the 
CCE shrink with increasing energy, the high momentum 
cusps resulting from inter-ellipse scattering move up with 
energy, and eventually merge. For the same reason, the 
low-momentum peaks resulting from intra-ellipse scat- 
tering disperse downwards with increasing energy. Since 
the negative energy CCE (not shown) expand with en- 
ergy, we find the opposite dispersion in this case as seen 
from Fig. 6(b). From the cuts along the antinodal direc- 
tions shown in Figs. 6(c,d), we see that the QPI response 
is actually strongest in this direction in agreement with 
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FIG. 6: QPI line cuts along the (110) (a,b) and (100) 
(c,d) directions for the pure AF case with M — t. Pan- 
els (a,c) [(b,d)] correspond to positive [negative] energies at 
\uj\/t = 0.01,0.075,0.15,0.225,0.3,0.375,0.45,0.5 (bottom to 
top). For clarity the curves are displaced by 0.15. 



the QPI maps in Fig. 5 

Finally we note that in samples with inhomogeneous 
coexisting regions of dSC and AF, one way to determine 
whether which contribution of the QPI signal originates 
mainly from the dSC or the AF regions, is to compare 
the result to the bias reversed QPI map. In the pure 
dSC state the same positions of q- spots should appear 
but with different weight, whereas in the AF phase the 
arcs in the QPI maps have also dispersed resulting in new 
locations of the cusps in e.g. the line cuts shown in Fig. 
6. Below, we study the simpler case of a homogeneous 
coexisting phase of dSC and AF long-range order. 



V. COEXISTING SUPERCONDUCTING AND 
ANTIFERROMAGNETIC PHASE 

In the homogeneous coexisting state (AF+dSC) with 
A, M 7^ 0, the d-wave gap collapses the arcs of spectral 
intensity at the Fermi level shown in Fig. 4(a,b) to nodal 
points. This robustness of the nodal points to antiferro- 
magnetism away from half-filling follows directly from the 
fact that (7r,7r) does not nest the nodal points. Some 
representative CCE for the coexisting phase are shown in 
Fig. 7. At low energies the CCE are again reminiscent of 
dispersing bananas, but now shifted off the normal state 
Fermi surface, and exist also in the shadow band outside 
the RBZ. This implies the existence of shadow banana 
tips and shadow QPI peaks, whose weight will however 
generally be suppressed due to coherence factors. 
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FIG. 7: (Color online) Contours of constant energy for the 
coexisting phase with A = 0.6t and M — t plotted at energies 
cj/t = 0.075, 0.15, 0.225, 0.3, 0.375, 0.45, 0.5, 0.6. The banana 
tips reach the AF zone boundary at a; = Ao = 0.33t when 
M = t. 
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FIG. 8: (Color online) QPI maps versus Qx/tt and Qy/Tr in the 
coexisting AF and dSC phase in the presence of single point- 
like nonmagnetic impurity. Energies uj are as shown in the 
Figure titles. This figure can be directly compared to Fig. 2. 
Note the disappearance of the spot-like feature corresponding 
to octet vector qr above the critical energy Ao = 0.33t. 



Once the banana tips reach the AF zone boundary, 
the CCE become similar to the pure AF case. In Fig. 
8 we show the results of the QPI patterns in the coex- 
isting phase. This figure can be directly compared to 
Fig. 2. We clearly see the mixing of the AF and dSC 
QPI features and the importance of the crossover scale 
set by the energy Aq where the CCE reach the AF zone 
boundary: at Cc; < Aq {u > Aq) the QPI is dominated by 
the pure dSC (AF) phase with the addition of possible 
shadow band features resulting from scattering involv- 
ing the bananas outside the RBZ. Peaks generated by 
shadow bands are particularly evident in the QPI image 
at Cc; = 0.225t in Fig. 8, and constitute a clear signature 
of the coexisting phase. For BSCCO, however, it seems 
likely that the disorder is simply too strong for these ad- 
ditional peaks to be observed [see next section]. In Fig. 
8 we have circled the position of the qr peak from the 
pure dSC phase [see also Fig. 2]. This peak appears to 
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FIG. 9: (Color onhne) QPI line cuts along the (110) (a,b) 
and (100) (c,d) directions for the pure coexisting AF and 
dSC phase with M — t and A = 0.6t. Panels a,c 
(b,d) correspond to the following positive (negative) energies: 
\uj\/t = 0.01, 0.075, 0.15, 0.225, 0.3, 0.375, 0.45, 0.5, 0.55, 0.6 
(bottom to top). The line cuts close to the crossover region 
ix) — Ao are indicated by red dashed lines. For clarity the 
curves are displaced by 0.15. 



be extinguished when crossing the AF zone boundary at 
cj = Aq. For the band parameter and superconducting 
gap A = 0.6t used in this paper, we have Aq = 0.45t for 
M 0.5t, and Aq 0.33t for M = t. 

Figure 9 shows the nodal and antinodal line cuts with 
the red dashed lines indicating the crossover region near 
Aq. By comparison to Figs. 3 and 6, the resemblance 
to the dSC and AF state is again clear for a; < Aq and 

> Ao, respectively. Note that despite the simplicity of 
the present model, it is evident from Figs. 8 and 9 that 
at a; > Aq the FT-STS maps are dominated by slowly- 
dispersing peaks in the antinodal direction qualitatively 
similar to the experimental data.^^ These maxima dis- 
perse, however, more than those in Ref. 36 and also ap- 
pear to be more arclike than the experimentally observed 
q^ and qg peaks. This may be because our simple model 
assumes a standard d-wave gap, whereas ARPES experi- 
ments reveal a gap which varies little near the antinodal 
points. We clearly also do not have in our model ad- 
ditional charge order leading to the nondispersive QPI 
peaks. 



VI. EFFECTS OF DISORDER 

So far the discussion has focussed on pure phases, dSC, 
AF, and AF+dSC in the presence of a single point-like 
nonmagnetic impurity. Clearly this situation is rather 
simplistic and should be extended to describe the QPI in 
a disordered short-range AF cluster glass phase as rele- 
vant for LSCO, and possibly also BSCCO. One way to 
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model short-range AF correlations is by assuming that 
the AF ordering vector Q follows a Lorentzian probabil- 
ity distribution pq peaked at (tt, tt) and with a broad- 
ening given by the inverse of the AF correlation 
length. The main effects of disordering with such a dis- 
tribution was studied in Ref. 67: only the part of the 
hole pockets which lie outside the RBZ is shifted. This 
means, for example, that scattering events involving the 
outer shadow band bananas [see Fig. 7] will be smeared 
out by this kind of disorder averaging, rendering the low- 
energy FT-STS patterns virtually identical to those of 
the pure dSC state. Nevertheless, modeling the disorder 
by a Pq distribution appears questionable for the CSG 
phase since it does not properly treat the spatial inhomo- 
geneity of the spin glass, and misses e.g. additional low- 
energy states existing at the boundary regions between 
magnetic and superconducting domains.''^ Therefore, it 
would be interesting to compare with quasiparticle inter- 
ference patterns from more realistic real-space disorder 
configurations similar to those produced in Refs. 51- 
54,56,61,62. The results presented in this paper should 
be helpful in providing an interpretation of results from 
future more realistic disorder modeling. 

VII. CONCLUSIONS 

We have studied quasiparticle interference phenomena 
in a dSC phase with coexisting short-range AF order, 
as found, e.g. in the cluster spin glass phase of under- 
doped cuprates. To do so, we have assumed that much of 
the qualitative physics may be captured by studying the 
analogous problem with long-range AF order. In partic- 
ular, we have calculated the QPI patterns arising from 
scattering from a single point-like impurity in the case 
of a pure metallic AF phase and a coexisting phase of 
both AF and dSC order. Due to different coherence 
factors in the AF and dSC phases, the QPI maps are 
dominated by peaks (arcs) in the dSC (AF) phase, re- 
spectively. In the case with coexisting order, low-energy 
quasiparticles propagate on contours which resemble the 
pure superconductor, so dispersive, localized interference 
spots similar to those predicted by the octet model for 
the pure dSC system are recovered. When the tips of 
these contours reach the AF zone face, however, the con- 
tours change abruptly to those characteristic of the pure 
AF, and the dispersing low-energy localized interference 
peaks are then extinguished, as reported by Kohsaka et 
al.^^ We expect that the remaining arc-like features at 
energies > Aq are suppressed by the short-range na- 
ture of the true magnetic order in these systems. 

The net result is a system where the low-energy quasi- 
particle interference features resemble those of the opti- 
mally doped superconducting samples, namely they dis- 
perse according to the octet model. For energies above 
the critical energy Aq, these localized spot-like features 
in momentum space effectively disappear, as seen in ex- 
periment. Although the formation of AF long-range or- 



der itself may be viewed as a coherent multiple scatter- 
ing process, we see from the study of the coexisting AF- 
dSC system that the quasiparticle states near the antin- 
ode are not destroyed by broadening in this process, but 
simply folded back. Thus in the realistic situation with 
short-range AF order, some scattering from AF modula- 
tions will occur, but still well-defined quasiparticle peaks 
should remain, as seen in ARPES. This resolves an appar- 
ent paradox in the comparison of the two experimental 
techniques in their view of the antinodal quasiparticles. 
It is the localized QPI q-space structures due to these 
states which are destroyed, not the states themselves. 

We emphasize that the formalism presented in this pa- 
per is essentially identical for any competing order sce- 
nario with ordering wavevector near (7r,7r). We believe 
that the bulk of the experimental evidence support the 
identification of the competing phase with short-range in- 
commensurate AF order or quasistatic fluctuations from 
incipient order, as observed in /iSR. We have therefore 
discussed the extinction of QPI above a critical energy in 
this context, but other explanations with a similar struc- 
ture may be possible. 

In addition, we note that realistic simulations of the 
QPI patterns require a more sophisticated description of 
disorder than that adopted here. It is thought that the 
potential produced by a single defect has components 
in at least the screened Coulomb and pairing channels, 
and some account of the latter is known to be necessary 
to reproduce the correlations between dopant position 
and gap size observed in the BSCCO system.-^'' We have 
neglected these details here, but they are discussed e.g. 
in Ref. 19. Nevertheless, the patterns produced by our 
single potential scatterer were found to be sufficient to 
describe the basic phenomenology of Kohsaka et al.^^ for 
the octet qy peak. It will be interesting to see if more 
sophisticated simulations can also reproduce the behavior 
of the other octet vectors and their weights correctly. 

One aspect of the current picture which remains un- 
clear is the extent to which static order is required. The 
QPI extinction phenomenon observed by Kohsaka et al.^^ 
is observed in BSCCO samples with doping levels 6-19%, 
i.e. including optimal doping where recent neutron ex- 
periments have claimed a 32 meV spin gap.*"^ Still, near 
(7r,7r), some intensity remains below this energy, and 
/iSR continues to indicate frozen magnetic order even 
in these samples at low temperatures.^^ It may be that 
the magnetism in this system is simply too disordered 
to be seen by current neutron experiments. A more so- 
phisticated treatment of short-range order is required to 
address these questions. 
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